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Abstract. Experimental data from the Relativistic Heavy Ion Collider (RHIC) 
suggests that the quark gluon plasma behaves almost like an ideal fluid. Due to 
its short lifetime, many QGP properties can only be inferred indirectly through a 
comparison of the final state measurements with transport model calculations. Among 
the current phenomena of interest are the interdependencies between two collective flow 
phenomena, elliptic and triangular flow. The former is mostly related to the initial 
geometry and collective expansion of the system whereas the latter is sensitive to the 
fluctuations of the initial state. For our investigation we use a hybrid transport model 
based on the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) transport 
approach using an ideal hydrodynamic expansion for the hot and dense stage. Using 
UrQMD initial conditions for an Au-Au collision, particles resulting from a collision 
are mapped into an energy density distribution that is evolved event-by-event with a 
hydrodynamic calculation. By averaging these distributions over different numbers of 
events, we have studied how the granularity/smoothness of the distribution affects the 
initial eccentricity, the initial triangularity, and the resulting flow components. The 
average elliptic flow in non central collisions is not sensitive to the granularity, while 
triangular flow is. The triangularity might thus provide a good measure of the amount 
of initial state fluctuations that is necessary to reproduce the experimental data. 
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1. Introduction 

Higher order flow coefficients have recently been recognized as a new observable to 
gain information about the creation of the quark gluon plasma in relativistic heavy ion 
reactions and its properties [U El [3j HJ [5] [6]. Elliptic flow, the second order Fourier 
coefficient of the azimuthal distribution of the particles in the final state, has been 
investigated in detail during the last decade, since it is the crucial observable that is 
responsible for the proof that the quark gluon plasma is a nearly perfect liquid [?]. 
Pressure gradients translate the initial state coordinate space eccentricity to the final 
state momentum space ellipticity and this connection is affected by the viscosity and 
the equation of state. 

Triangular flow and even higher moments of the Fourier decomposition have been 
measured during the last year by various experimental collaborations at RHIC and LHC 
[HJ |9j [10] . The higher odd anisotropic flow coefficients require the treatment of event- 
by-event fluctuations, since they vanish averaged over events. A lot of recent theoretical 
development has been aimed towards a better understanding of the correlations between 
initial and final state. Most of these models work in a simplified framework, using e.g. 
linearized hydrodynamics or rotated averaged initial conditions [UJ [12], [13] [H] . 

The initial state of a heavy ion reaction that is needed to run a (viscous) 
hydrodynamic evolution is usually modeled following a Glauber or CGC approach |15j . 
Also more realistic initial conditions have been investigated using event generators like 
NEXUS, URQMD, etc PS] [17]. Some people have also studied simplified models like the 
one-tube model [18] or just a statistical number of hot spots to get a more systematic 
handle on understanding initial state structures [T9] [20]. 

Between all these different approaches there seems to be agreement, that the flow 
coefficients can be used to constrain the initial state dynamics and have the potential to 
sort out different models [21J, [22]. The ultimate goal is to understand the initial energy 
deposition which has fundamentally to do with the distributions of the nucleons/partons 
in the incoming nuclei and the interactions they are undergoing. 

The systematic study that is presented in this paper constitutes one step on the 
way to constrain the initial state granularity. By using a hybrid transport approach that 
is based on the Ultra-relativistic Quantum Molecular Dynamics including an (3+1) 
dimensional ideal hydrodynamic expansion, we demonstrate that triangular flow is 
directly related to the amount of fluctuations in the initial state. By using a full event- 
by-event setup we can calculate the triangular flow via the event plane method as it is 
measured by the PHENIX collaboration and put some first constraints on the initial 
state granularity. 

In Section [2] the hybrid model is described and it is explained how different 
granularities are obtained by initially averaging over different numbers of events. The 
following Section [3] demonstrates that the chosen setup does not change the bulk 
properties of the system by looking at pion and kaon transverse mass spectra and elliptic 
flow in non-central collisions. After that, we come back (in Section H]) to the initial state 
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eccentricity and triangularity distributions as a function of granularity leading directly 
to the resulting anisotropic flow coefficients in Section [51 In the end, we will draw 
conclusions from the presented results and outline further research. 

2. Initial State Granularity 

The initial state profile and its granularity are mainly influenced by the following three 
things: the shape of the incoming nuclei, the interaction mechanism and the distribution 
of binary collisions and the degree of thermalization. In our approach all of these 
are given by a dynamic transport model, namely Ultra-relativistic Quantum Molecular 
Dynamics (UrQMD) [23j EH EH]- The incoming nucleons are sampled according to 
Woods-Saxon profiles and the initial nucleon-nucleon scatterings and non-equilibrium 
dynamics proceed according to the Boltzmann equation. After the two nuclei have 
passed through each other, at t start = 0.5 fm, local thermal equilibrium is assumed to 
perform the transition to the ideal hydrodynamic description. 

The particle degrees of freedom are represented by three-dimensional Gaussian 
distributions in the following way 



to obtain energy, momentum and net baryon density distributions that are smooth 
enough for the hydrodynamic evolution. Here, e is the energy density at position (x, y, z) 
that a particle with energy E p at position (x p ,y p ,z p ) contributes. The Gaussians are 
Lorentz contracted in z-direction by 7 2 to account for the large longitudinal velocities. 
Only the matter at midrapidity \y\ < 2 is assumed to be locally equilibrated, whereas 
the other hadrons are treated in the hadronic cascade. 

A similar procedure is often employed in Monte Carlo approaches for hydrodynamic 
initial conditions. The width of the Gaussians representing the individual particles or 
tubes that lead to hot spots in the initial state is then varied and the granularity changes. 
The disadvantage of adjusting this parameter is that the overall features of the initial 
state change as well, e.g. the entropy density or the maximum value of the profile and 
need to be re- adjusted. In [26] we have shown how bulk observables are affected by 
variations of the obvious initial state parameters, the starting time and the Gaussian 
width. 

For the present systematic study we have chosen a different path: By averaging over 
different numbers (n — 1 - n — 25) of initial events at the same beam energy, impact 
parameter and collision system, the bulk properties of the system are kept constant, 
while the granularity is varied. The disadvantage of the averaging procedure is that the 
direct connection to the initial state physics is lost. Still, this setup is a powerful tool 
to show systematic dependencies on the amount of initial state fluctuations. 

Fig. [1] and El show examples for initial states that are generated by averaging over 
1, 5 or 25 events. The maximum energy density in the transverse plane in central and 
non-central collisions is similar in the three cases and the granularity differs from a 
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Figure 1. Energy density distributions in the transverse plane for central (6 = 2 
fm) Au+Au collisions at ^/snn = 200 GeV with three different granularities. The left 
picture depicts the profile from one event, the middle one has been averaged over 5 
events and the right panel is a smooth average of 25 events. 



maximum amount of fluctuations using the full event-by-event setup (n = 1) to pretty 
smooth initial conditions (n = 25). We have generated 300 of each of these initial 
profiles for 6 different granularities for Au+Au collisions at a/snn = 200^4 GeV at two 
impact parameters. 




Figure 2. Energy density distributions in the transverse plane for mid-central (6 = 7 
fm) Au+Au collisions at ^/snn = 200 GeV with three different granularities. The left 
picture depicts the profile from one event, the middle one has been averaged over 5 
events and the right panel is a smooth average of 25 events. 

These individual initial conditions are propagated through a full (3+1) dimensional 
ideal hydrodynamic evolution [271 128] and a subsequent hadronic afterburner. For the 
equation of state and the freeze-out criterion the most successful values found in [2H] 
have been used (DE-EoS and e <= 1.022 GeV/fm 3 ). 

The event-by-event calculation provides the full final phase-space distribution of the 
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hadrons that are also measured in experiments and therefore we are able to employ the 
same event plane measurement method as the PHENIX collaboration has used for their 
V3 analysis. This allows us to directly compare to experimental data and investigate the 
sensitivity to the granularity. 



3. Spectra and Elliptic Flow 

In Fig. [3] a comparison between the hybrid approach for different granularities and the 
experimental data for basic bulk observables is shown. The left hand figure depicts 
transverse mass spectra for pions and kaons in central collisions and the right hand 
figure displays elliptic flow for charged particles as a function of transverse momentum 
in mid-central collisions. 
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Figure 3. Transverse mass spectra for ir~ and K + (left) and elliptic flow as a function 
of transverse momentum for charged particles (right) at midrapidity (|y| < 0.5) in 
central/mid-central (6 = 2 fm/6 = 7 fm) Au+Au collisions at ^/snn = 200 GeV 
from the hybrid approach for different granularities compared to experimental data 

[SSlEolEUIH]. 



The elliptic flow has been calculated using the event-plane method including 
corrections to remove auto-correlations and resolution as described in [I] which matches 
the PHENIX analysis [8]. The major observation is that the transverse mass spectra 
and the elliptic flow do not depend on the granularity. This result proves that our 
setup leaves the bulk properties unchanged [32] and is in reasonable agreement with the 
experimental data. 



4. Eccentricity and Triangularity Distributions 

To exploit the potential sensitivity of final state flow coefficients on the granularity 
let us first have a look at the initial state coordinate space Fourier coefficients. The 
eccentricity and the triangularity have been calculated from the center transverse slice 



Initial State Granularity 6 

(z = Ofm) of the initial energy density distribution using the following formula: 
J (r n cos(n0)) 2 + (r n sin(n0)) 2 

£ " = H • (2) 

(, ) depicts an average over the cells of the hydrodynamic grid weighted by the energy 
density of the respective cell and r and <fi are polar coordinates. 
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Figure 4. Eccentricity distributions for different granularities in central/mid-central 
(left: b = 2 fm and right: b = 7 fm) Au+Au collisions at ^/snn = 200 GeV. 

In Fig. H] the distribution of the eccentricity for 6 different granularities in central 
(left) and mid-central (right) collisions is shown. In the later case, the mean value of 62 
is remarkably stable and the distribution gets only more and more peaked for smoother 
initial conditions. This behavior is consistent with the constant elliptic flow result in 
the previous Section. For central collisions the major source of elliptic flow are initial 
state fluctuations and not the geometry given by the reaction plane as in non-central 
collisions, therefore the mean value of eccentricity decreases for smoother more central 
collisions. 

Triangular flow and the initial state triangularity are only generated by initial 
state fluctuations independent of centrality (as long as the collision is not ultra- 
peripheral [33]). Therefore, the mean values and the general behavior of the triangularity 
distributions shown in Fig. [5] are only slightly dependent on centrality. The triangularity 
decreases to almost zero for smooth initial conditions and shows a clear sensitivity to 
the granularity. These observations can be summarized in the following statement: The 
mean value of the triangularity and the fluctuations of the eccentricity are sensitive to 
the initial state granularity. Since fluctuations of elliptic flow are harder to measure 
and to compare consistently to model calculations, it seems more promising to look at 
triangular flow to constrain the amount of fluctuations in the initial state. 
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Figure 5. Triangularity distributions for different granularities central/mid-central 
(left: b = 2 fm and right: b = 7 fm) Au+Au collisions at ^/snn = 200 GeV. 



5. Anisotropic Flow Results 

In this Section we present results for the averaged t>2,3 coefficients for charged particles 
calculated for different granularities (n = 1,...,25) via the event plane method. The 
results in Fig. [6] have not been corrected by the resolution factor since the resolution 
gets very small for smooth initial conditions (with the exception of i>2 in non-central 
collisions), where the flow coefficients almost vanish anyhow. 
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Figure 6. Anisotropic flow for different granularities in central/mid-central (left: 
b = 2 fm and right: 6 = 7 fm) Au+Au collisions at ^/snn = 200 GeV. 



The clear dependency of the mean values of triangular flow and elliptic flow in 
central collisions on the number of events over which the initial average has been 
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performed proves that one can constrain the granularity of the initial state by comparing 
flow coefficients of calculations to experimental data. 




Figure 7. Triangular flow versus transverse momentum for different granularities 
central/mid-central (left: b = 2 fm and right: b — 7 fm) Au+Au collisions at 
Tsi^ = 200 GeV compared to PHENIX data [8]. 



To get a first impression on the amount of fluctuations that is consistent with the 
PHENIX data on triangular flow [8] , the transverse momentum dependence of triangular 
flow of charged particles has been calculated in the hybrid approach (see Fig. [7]). Up 
to a value of n = 5 the resolution of the event plane has reasonable values and the 
calculation has been corrected by the resolution factor. For the two smoother cases 
again the resolution gets too small and the triangular flow is practically zero. This first 
comparison to data indicates that the full event-by-event setup is close to the amount 
of initial state fluctuations that is necessary. 

Since additional viscosity during the hydrodynamic evolution would dilute 
fluctuations faster this is a lower bound on the initial state granularity. If there would 
be a similar event-by-event approach available with a viscous hydrodynamic evolution 
one could first fit the viscosity to elliptic flow in non-central collisions and then tune the 
initial state granularity to match the triangular flow data. 



6. Conclusions and Outlook 



In this systematic study of we have shown that the mean value of triangular flow can 
be used to constrain the initial state granularity in heavy ion reactions. A (3+1) 
dimensional hybrid transport approach that is based on UrQMD including an ideal 
hydrodynamic expansion has been employed. The granularity of the initial state can be 
varied by averaging over different numbers of events before generating the energy density 
distribution for the hydrodynamic evolution. This setup has the main advantage that 
the bulk properties of the system are not affected and match the experimental data. It 
is shown that the full event-by-event setup gives currently the best description of the 
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available PHENIX data. 

One future direction of research is to define a measure of the initial state granularity 
with a more intuitive interpretation, that can be applied to arbitrary initial state profiles. 
The wealth of experimental data on higher order flow coefficients in turn will help to 
pin down this granularity measure and define the initial state profiles. Our plan is to 
perform a multi-paramater sensitivity analysis using sophisticated statistical methods 
for that purpose [53] . 
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